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O Abstract 

, In this article we study a problem related to the first passage and inverse first 

passage time problems for Brownian motions originally formulated by | Jackson,] 



> 
in 



Kreinin, and Zhang (2009)] Specifically, define tx = mf{t > : Wt + X < b{t)} 
where Wt is a standard Brownian motion, then given a boundary function 
b : [0, oo) M and a target measure fi on [0, oo), we seek the random variable 
X such that the law of tx is given by ^. We characterize the solutions, prove 
uniqueness and existence and provide several key examples associated with the 
linear boundary. 

1. Introduction 

In this paper we examine a problem related to both the first-passage and inverse passage time 



problems. It was originally formulated by Jackson, Kreinin, and Zhang (2009) for the case of 



g a linear boundary. In general, let (fi.P.F) denote a eomplete probabiUty spaoe, and let F = 

{Ft}o<t<T denote the natural filtration generated by the standard Brownian motion Wt- Consider 
^ a deterministic function b : [0, oo) M with 6(0) > — oo (this will represent the stopping boundary 

^ for the Brownian motion) and a random variable X > 6(0) (this will represent the randomization of 

the initial starting point of the Brownian motion) . We will assume that X spans the cj-algebra Fq 
and is independent of the Brownian motion. Define the stopping time of the randomized Brownian 
motion hitting the boundary 6(t) as follows 

TX = inf{ t > ; Wt + X < b{t) } . (1) 

Without loss of generality we can take 6(0) = so that X is non-negative (see Figure [T|. Given the 
above, we are interested in the following problem: 
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Figure 1: A sample path of the randomized Brownian motion W{t) +X. g(x) is the density of the starting point and 
Tx is the first passage time of the randomized Brownian motion to the boundary b{t). 

Definition 1.1. (Randomized First Passage Time Problem (RFPT).) Given a boundary function 
h : [0, oo) M, and a probability measure /x on [0, oo), find a random variable X such that /j, is the 
law of the randomized first passage time tx ■ 

Although the RFPT problem is stated in full generality, we will assume that b{t) is continuously 
differentiable in order to ensure that the conditional distribution of t\x=x has a continuous density 



function f{t\x) (see Peskir (2002)). Furthermore, we will assume that fi is absolutely continuous 



with density function f{t) which we refer to as the target density. Then the RFPT probelm is 
equivalent to solving the Fredholm integral equation: 



fit\x)gix)dx = fit) . 



(2) 



Here, the conditional density of the hitting time f{t\x) acts as the kernel in the integral equation. 
In general we seek solutions to ^ in the class of functions Q := { g ; J^_|_ \g\ < oo or \g\ < K }; 
however, our main focus is on the class of density functions on the positive real line which is included 
in Q. 

In this paper we address the questions of existence (see Theorem [T]) and uniqueness (see Propo- 
2.1 and 2.2 ) of the distribution of the random initial point X. Theorem [T]is one of the main 



sitions 



results of this work and it implies the existence of a random variable X with moment generating 
function g given by 



9(a) 



^-ab{t)-a'^t/2 



f{t)dt 



(3) 
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This existence result is very powerful, and provides a relationship between the Laplace transform 
of the initial point and a boundary specific integral transform of the target density. For particular 
boundaries, this integral transform is of a standard type. For example, if the boundary is linear 
h{t) = /xt, then the integral transform in the rhs of ([s]) is simply the Laplace transform at the point 
a// + \oP'- However, in general the integral transform may not be easily identifiable; nonetheless, 
one can in principle compute the integral transform numerically and then numerically invert g to 
obtain the density. 

Note that, since the boundary h{t) defines the conditional distribution of the stopping time f{t\x) 
uniquely, when X is a discrete random variable the distribution of the unconditional stopping time 
will be a mixture of the conditional densities as seen from ([2|. The reverse statement, however, 
does not hold. We show this in the case h{t) = fit. Given the well known fact that continuous 
distributions on the positive real line can be arbitrarily well approximated (at least point-wise) by a 



mixture of Gamma distributions (see e.g. Tjims (1995) ), we derive analytical formula for the density 



function of X when the target distribution of tx is a finite mixture of Gamma distributions. This 
formula allows us to extend our analytical results for g to a very large class of target distributions. 

When X is non-random, the RFPT problem reduces to the classical first-passage time (FPT) 
problem for Brownian motion if we take the boundary function 6(.) as input and seek the density 
function /(.). The FPT problem has a long history and dates back to the work of A.N. Kolmogorov, 
A. Khinchin, L Petrovsky and P. Levy. A good summary of the their early work can be found 
in Khinchine (1933) The available closed form results appear to be sparse, fragmentary and 



essentially confined to the linear, quadratic (see Salminen (1988) I and square root (see Novikov 



(1981)) boundaries. A unified approach for the derivation of the analytical results for these three 



cases was presented in Jaimungal, Kreinin, and Valov (2009) based on the Fredholm equation 



roo 

-ab{t)-ta'^ /2 . 



r e-"''^'^-'°'/V(t)dt = 1 , (4) 
Jo 



where / is the density function of the first passage time of the Brownian motion to the boundary 
b{t). Equation ^ holds for all a and continuous functions 6(.), 6(0) > — cxd, satisfying 

lim(6(t) + at) > -oo . (5) 

oo 

This condition simply states that the drifted boundary is uniformly bounded below. 



Jaimungal, Kreinin, and Valov (2009) also derive a class of Volterra Integral equations of first 
kind which generalizes and unifies all previous known FPT/IFPT integral equations. A key element 
in that work was the construction of the following new class of martingales. 
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Proposition 1.1. The process 

Zs := m{s,Ws;p,t), m{s,w;p,t) = _ g)(p+i)/2 ' 

is a real valued martingale on s € [0,t) for allp,y E M, t > 0. Here Dp{z) is the parabolic cylinder 
function - for some of its properties see Appendix \B.l 



Armed with this martingale, the authors then derive the class of Volterra integral equations 

m{t,0;p,t)= m{s,b{sy,p,t)f{s) ds , (7) 
Jo 

for y < b(t), by invoking the optional sampling theorem on the martingale Zs- Furthermore, by 
passing to the limit y | b{t), the authors show that the resulting class of Volterra equations contains 



the class of equations of Peskir (2002) and show that the FPT density is the unique continuous 
solution to any member of this class under certain regularity conditions on the boundary. 

When X is non-random and we take the density function /(.) as an input, the RFPT problem 
reduces to the inverse first-passage time (IFPT) problem which seeks the boundary function 
The problem was first posed by A. Shiryaev in 1976 for the case of the exponential probability 



density. An early paper by Anulova (1980) deals with the existence of some stopping times for a 
given distribution, however, these stopping times are not of the form ([T| for some function b. One of 
the main contributions to the description of the inverse problem thus far, was provided by [Chadam 



Chen, Cheng, and Saunders (2006) where they demonstrate the existence of a unique viscosity 



solution for the IFPT problem, from a PDF perspective, and describe the small time behaviour of 
the boundary function. Note that in the context of the IFPT problem the integral equations ([7| are 
nonlinear. 

Randomizing the starting point of the Brownian motion allows us to bypass both the first passage 
time and inverse first passage time problems. This is achieved by assuming the pair (b, f) as given 
while trying to match the density / by randomizing with a density function g. Furthermore the 
randomization allows the distribution of the hitting time tx, and thus of the stopped process WV^j 
to probe a much wider class of distributions than if the Brownian motion starts at a fixed point. 
For example, in the linear boundary case, with slope equal to one, we have the relation Wr^ = 
Tx — X. If we can imply the distribution of X from any distribution of Tx, then we may be able 
to describe the resulting class of distributions for WVx • This is tantalizingly close to the statement 



of Skorohod's embedding problem (see Skorohod (1965) I, where one seeks a stopping time r such 
that the stopped Brownian motion Wr has a given distribution. Here, our class of stopping times 
Tx is generated by fixing the boundary b and randomizing the starting point X, and through the 
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connection Wr^ = bijx) the distributions of the stopped Brownian motion and the stopping time 
is clear. 

If we attack the RFPT problem by seeking a direct solution to ^ it may seem, at a first glance, 
that this could be a formidable task since the kernel of this Fredholm equation, f{t\x)^ is unknown 
for most boundary functions. However, as we will discover below, the problem is simpler than 
both the FPT and IFPT problems and we obtain analytical and semi-analytical results for certain 
transforms of the matching distribution using the Volterra and Fredholm integral equations derived 
J aimung Kreinin, and Valov (2009) [ 



m 



The remainder of the paper is organized as follows. In Section 1, we assume the existence of a 
function g which solves ^ and, under certain conditions on the boundary 6, derive unique integral 
transforms of g. The integral transforms lead to uniqueness of the solution and provide us with a 
way to compute g analytically. In Section 2, we address the question of existence of the random 
variable X. In Section 3, we look at the linear boundary case and compute g analytically for a class 
of target distributions. Moreover, we present a number of examples for the pair {g,f). In Section 
4, we examine the effect of an affine boundary transformation on the target density. For the case 
of the linear boundary, we explore the relationship between the target densities corresponding to 
different slopes assuming the same density function g. Furthermore, we analyze the implications of 
a certain scaling property of the boundary and in the linear boundary case we look at an extension 
of the RFPT problem when both the intercept and the slope are random. Finally, we end with some 
concluding remarks in Section 5. A number of technical proofs are delegated to the appendix. 

2. Uniqueness 

The integral equations ([7| take on a particularly useful form when p is restricted to the integers. 
In this case, the parabolic cylinder function Dn is related to the Hermite polynomials Hn as follows 
(n e Z): 

Dn{z) = (z/V2) . (8) 

This allows equations ([7|, for p = n, to be rewritten as: 



for any y < h{t) — x. Here we have the equation in the RFPT form by using the conditional 
hitting density f{t\x). Note that since the Hermite polynomials form a complete orthogonal basis 
in L^(M,e~^^) with respect to the standard normal distribution, ([9| allows for a unique series 
representation of the density function g whenever it exists. Here we seek uniqueness for the larger 



6 S. Jaimungal, A. Kreinin, and A. Valov 

class Q which contains, but is not Hmited to, the class of densities. This leads us to the following 
results. 

Proposition 2.1. Suppose b : [0, oo) i— > M is a continuous function and there exists a t > such 
that b{t) > 0. Then, if Q) has a solution g G G P| L^(M, e~^^), it is unique and it is given by 

a(-) = E -^H^r.{x/V2i) , (10) 

^ n=0 

for any t > such that b{t) > 0, and where 
Proof. See Appendix A □ 

This solution need not be a true density function. However, in the case when g' is a density 
function then it is the unique solution to the randomized FPT problem. Interesting, the above 
series representation holds for all t > such that b{t) > 0, yet the solution g(x) is independent of 
the specific choice of t. 

Next we examine the Laplace transform of g using the Fredholm equation of the first kind 
Using ([2]) we obtain the following result. 

Proposition 2.2. Suppose b : [0,oo) 1-^ M is continuous and satisfies condition (jsj). Then, if ^ 
has a solution g £ G, it is unique and its Laplace transform is 

~g{a) = / e-^'^'^-'^'/^f{t)dt . (12) 
Jo 

Proof. Assume that equation ^ has a solution g £ Q. Under the condition (|5|, equation Q 
holds for all a > and (conditional on X = x) we have 







Multiply both sides by g{x) and integrate out x. By Fubini's theorem, the order of integration can 
be exchanged since 

/ / \g(x)\e-''^^''^-*°'"l'^f{t\x)dtdx = / e-"n5(a;)|dx < 00 . 

JO JO JO 

Finally, using (|2]), we obtain the Laplace transform of g given in (12). Uniqueness then follows from 
the uniqueness of the Laplace transform. □ 

As a demonstration of the applicability of the above result, suppose b{t) = ^/t and let / be the 



unconditional density. Then (12) becomes 



~9{a) = / e-^^'-'^'/y{t)dt , 
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for all a > 0. Multiplying both sides of this equation by a^"^, p > 0, and integrating out a we 
obtain 



/•oo 

Jo 



oo />oo 

Jo 

oo POO POO 

t~P/^f{t) / uP-ie-"-"'/2 ^^^^ ^ ei/^r(p)D_p(l) / t-P^^f{t) dt 
Jo Jo 

= eV4r(p)D_,(l)7(l-p/2) . 

Here, the substitution u = a\/t was used in the second equality, and / is the Mellin transform of /. 
Thus the Mellin transform of the Laplace transform of g, denoted g is given by 

f(p) = ei/^r(p)D_p(l)/(l-p/2) , 

provided that /(I — p/2) exists for a non-empty set of positive real values of p. 

3. Existence 



We saw from Propositions 2.1 and 2.2 that if the boundary is well behaved and there exists a 



solution g G ^ of ([2| then it is unique. Thus existence of a solution to the RFPT problem implies 



uniqueness if the boundary satisfies the hypothesis of Proposition 2.2 since any density function g 
belongs to the class Q. The question of showing the existence of a unique matching density g reduces 
to finding conditions under which Q has a solution. Given b{t), there may not exist a density g for 
every density function / satisfying d2|. An illustrative counterexample of existence is provided by 



the choice b{t) = and f{t) = Ae"^*, A > 0. Then ([12]) reduces to g{a) = f{a^/2) = 2A/(2A + a^). 



Consequently, g is the Laplace transform of g{x] A) = v2Asin(xv2A). Thus clearly, g & G but it is 
not a probability density function. 

A sufficient requirement for the existence of a density solution to ^ can be constructed based 
on Picard's Criterion (see e.g. Polyanin and Manzhirov (2008)] p. 578-583). However, such a 



construction is difficult since the kernel function f{t\x) is not, in general, known explicitly. The 
only notable exceptions being the linear, square-root and quadratic boundaries mentioned earlier. 
Furthermore we cannot guarantee that the solution is a density function. Consequently, we approach 
the question of existence in a probabilistic manner. To this end, we seek a random variable X such 
that 



Ei{f{t\X)) = f{t) , 



(13) 
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where the expectation is taken under a measure q with support on the positive real hne. In Hght 



of the Laplace transform given in (12), the question of existence is reduced to examining sufficient 
conditions under which this transform is a moment generating function of some random variable 
X. It is well known that moment generating functions are completely monoton^(see Feller (1971) |. 



Thus, if the function r : [0,oo) ^ (0,oo), defined as 

/>oo 

r(a) := / e-"^W~"'*/2_^(t)dt , (14) 







is completely monotone, then it is our candidate for a moment generating function. The following 
Lemma provides an alternative check point for checking complete monotonicity of r{a) and proves 
to be a useful tool. 

Lemma 3.1. Suppose b : [0, oo) h^M. is continuous, satisfies condition and 

e'^'-^f{t)dt < oo, V/3 > . (15) 
Then r(a) is completely monotone if and only if the pair (b, f) satisfies 







for all n > 1 and a > 0. 
Proof. See Appendix A □ 

Using the complete monotonicity property of r(a), we can now derive sufficient conditions for 
the existence of X. The following Theorem describes the properties of the functions b and / which 
guarantee the existence of X. 

Theorem 1. Suppose that 

1. b{t) : [0, oo) 1-^ M is continuous and satisfies condition ^ 

2. r{a) is completely monotone 

3. If 

POO 

/ e-"''W-"'*/2^(t)dt = (17) 
Jo 

holds for all q > 0, then z(t) is identically zero 

Then, there exists a random variable X with m.g.f. given by r{a) such that tx has probability density 
function f. 



Recall that a completely monotone function r has derivatives of all orders which satisfy {—l)"-£-i^r{a) > for all a > and 
all non-neeative integers n > 0. 
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Proof. Since r is completely monotone and r(0) = 1, by Bernstein's theorem (see Feller (1971) 
pp. 439), there exists a probability measure on [0,oo) with cumulative distribution function q such 
that 

r[a) = / e-^"dg(j;) . 
On the other hand since equation (|4| holds for all a > then we have 



oo 

e~ 





for all a > 0. Taking integrals on both sides of the above equation with respect to the function q 
and using Fubini's theorem we obtain 

r{a) = / e-"^dg(x) = / / e-''''^'^-''"'/'^ f{t\x)dtdq{x) 
Jo Jo Jo 

rca roo 

^-ab(t)-^H/2 / f{t\x)dq{x)dt . 

Jo 



The above relation implies that 



oo / poo 

g-«M^)-^t/2^ f{t\x)dq{x)-f{t)]dt = ^. (18) 



Since (|18|) holds for all a > 0, assumption 3 implies that 

f{t\x)dq{x) = fit) . 







Furthermore, integrating the above w.r.t. t on [0, oo) we have that Jq°° dq{x) = 1 since / is a proper 
density function. Therefore q defines a proper distribution function. Consequently, the equality 
E'^(/(t|X)) = /(t), where the distribution of X is given by g, holds and there exists a solution to 



the RFPT problem (13). □ 



The hypothesis of Theorem [T] imply the assumptions of Proposition 2.2 and therefore existence 



implies uniqueness. The assumption that (17) has only the trivial solution is certainly satisfied in 



the case b{t) = /it, fi > because, in this case, (17) implies that the Laplace transform of z{t) is 



zero and thus z{t) is identically zero. The assumption is also satisfied when b{t) = n\/t. In this case. 



applying the Mellin transform on both sides of (17), implies that the Mellin transform of z is zero 



and thus z is also zero. These two examples demonstrate that the class of boundaries for which (17) 
has only the trivial solution is non-empty. The complete monotonicity of the function r is harder to 
check in general. However, the case b{t) = fit and f{t) is the Gamma(a, b) density can be checked 
directly using the conditions of Lemma |3.1[ 



Corollary 3.1. Let b{t) = fit, fi > 0, and f{t) = with a > 2//x^. Then there exists a 

unique random variable X such that tx has the probability density f . 
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Figure 2: The probabilistic interpretation of equation (201. The Brownian path must first hit the boundary h before 
hitting zero. 



Proof. Clearly, b{t) satisfies the hypothesis of Proposition 2.2 which shows the uniqueness part. 
For proof of existence see Appendix A. □ 



Corollary 3.1 motivates the next section where we take a closer look at the linear boundary case. 

We end this section by showing that if there exists a unique solution g to the RFPT problem 
for a given boundary b{t) > 0, Vt > 0, and target density fb{t), then g is also a solution for the 
zero boundary with a particular target density /q. To show this relationship, consider equation 
When n = 1 (and y = —x) and using Hi{x) = 2x, (|9| reduces to 

fis\x)ds = I— ^.^ . (19) 



e 2(t-s)5(s) 



/o V2^(t - s)3/2 ^ ' ' V2^t3/2 
The right side can be recognized as the probability density of the FPT of x+Wt to the zero boundary. 
This observation admits a simple probabilistic interpretation of ( |19[ ): the process x + Wt first hits b 
at time Tb before it hits the zero boundary at time tq - see Figure |2j Thus, if there exists a g £ Q 
which solves ^ for the boundary b and unconditional density fb, then 



e 2(t-=)5(s) 



v27r(t-s)3/2 



Ms)ds = Mt) 



(20) 



where fb{s) = f{s\x)g{x)dx and fo{t) = /q 



g{x)dx. If b{t) > for ah t > then ( 19 ) and 



(20) hold for all t > 0. In fact fo{t) is a proper density function, as can be seen by integrating the left 



side of (20). Therefore, g G GG is furthermore a density function, the corresponding unconditional 



distributions of the first passage times to b{s) and to are related as in (20). Moreover, (20) implies 
that for every distribution fb, for which there is a matching distribution g, there exists a distribution 
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/o (given by the integral in (20)) such that the pair (/o,& = 0) has the same matching distribution 
g. A.S a result, the class of unconditional densities for the boundary b{t) = for which there exists 
a matching distribution is at least as large as the corresponding class of unconditional distributions 
for any boundary h{t) > 0. 

4. Linear Boundary 

We now focus on the case when the boundary is linear and seek explicit solutions to the RFPT. 
When the starting position X is non-random then the first passage time distribution of the Brownian 
motion to the linear boundary bt — X is well known to be inverse Gaussian and is explicilty 



X r_(6t-X^2 



/.|^(t) = -^exp^- ^"^ \ . (21) 



Jackson, Kreinin, and Zhang (2009) use this explicit form to demonstrate that the hitting time of a 



drifted Brownian motion with a random starting point can be Gamma distributed. In this section 



we corroborate this result based on our integral equation (12 ) and extend it to a class of distributions 
which are infinite linear combinations of Gamma distributions. 
Letting b{t) = fit, n > and denoting 

Tx,f, = mf{t > 0;Wt < nt - x} , 



then, for a > —fi (12) reduces to 



oo 



5(a) = / e-*^°'^+" ^'^f{t)dt = fiafi + a'/2) , (22) 
Jo 

where / is the Laplace transform of /, the distribution of the randomized stopping time tx,^- When 
/ is the density of the Gamma distribution then we can factorize f{afi + a'^/2) and write it as a 
product of two Laplace transforms of Gamma densities. As a consequence, 5 is a convolution of 
Gamma distributions. The same argument applies when / is a mixture of Gamma distributions. 
More formally, define the sequence such that inf„>i > 2//x^ and the sequence 6„ > 0, Vn G N. 
Define the class of densities, 

!oo 00 ^ 

/ : fit) = J]pn/a„,b„(t), Y.P^ = ^,Pn> 0, inf a„ > 2/^2 I ^ (23) 
n=l n=l J 

where fa„,b„ are densities of Gamma distributions with scale parameters a„ and shape parameters 
bn- Then we have the following result: 

Theorem 4.1. Let b(t) = ^t, /i > 0. Suppose tx,^ ^ f £ C . Then, the matching density g{x) is 
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given by 

g{x) = Y.P^ V(h\ fTT / 2 ^1 K~l/2{xy/JF^2/^n) ■ (24) 



Proof. First we show that (24) holds for a finite mixture of Gamma densitites. Let f{t) = 
Y2n=i Pnfa^,br,{t) , J2n=iPn = ^7 Pn > 0, where /a„,fe„ are Gamma densities with scale parameter On 



and shape parameter bn- Then, from (22), we have 

N 



g{a) = fifia + o? 12) = + a„(^a + jl))'^- 

n=\ 

N 

= Y,Pn{l + c+a)-^"(l + c-a)-'" (25) 

n=l 

where = ^a„(/iib y^^^ — 2/a„). For 3ft(c^) > and 9(c^) = 0, i.e. to be positive real numbers, 
we require that a„ > 2//i^, n = 1, A^. This restriction is already enforced by having f £ C. From 
(25), g{x) is a mixture of convolutions of Gamma r.v.'s and in particular 

9{x) = Vpn / f^+ hSu)f^- hJx-u)du 
I Jo 

N I / \ fe-.-l/2 

^ \fiT{e~^^ ( X \ 

where / is the modified Bessel function of the first kind. 

For an infinite mixture of gamma distributions the result now follows easily. For tx/j ~ / G C, 



substitute / in (12) and using Fubini's theorem we can exchange the integration and summation 
(since all quantities are positive) to obtain the above result. The condition inf„>i a„ > 2//x^ ensures 
that each Laplace transform in the infinite mixture is factorizable with real valued roots. □ 

Next we look at several simple examples for a finite and infinite Gamma density mixture. 

• Example 1: For = 1 and ai = 2//i^ we have cf = and 

~g{a) = (1 + {l/fi){fia + = (1 + a/fi)-^''^ . 

Thus g is the density of a Gamma(l//i, 26i) distribution. 

• Example 2: For = 1, ai = 2//u^ and bi = 1/2 then g{x) = ^e~^^, the density of an 
exponentially distributed random variable with ^. 

• Example 3: For bn = ^ and fln ^ 2//x^, n < N, using the equality Ii/2{u) = 2 sinh(n) / \/27rn, 
we obtain 



- 2 /an 
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• Example 4: For /x > 1, take oi = 2, bi = k/2 so that rx,^ ~ X^(^) a-iid 5 is given by 



\k/2) \ 2^/^F^ 



While direct inversion of (|22| could be complicated for a general density /, Theorem 4.1 provides 
a procedural approach for computing densities from the class C and their corresponding matching 



densities given by (24) by simply choosing the sequences a„, 6„, p„. Due to the restriction on the 
scale parameters (as in the class C) we could set a fixed scale parameter = 1/c, c < /U^/2 and 
choose a set of shape parameters 6„ and weights pn to match a particular density. In such cases, the 
unconditional density function / becomes 

oo 

/(t)=ce-*^^(ct)''-.-^ (26) 

This class of densities includes the non-central x^{m) distribution by choosing 6„ = m/2 + n, c = 
1/2, Pn = (n-ty. — where 5 is the non-central parameter. Some more general examples of 



distributions of the form (26) are given below. 

• Example 5: pn = e~° ' = v + n, an = '^/c. Then / and g are given by: 



fc=0 



k\T{v + k + l) ' 



v+l/2 



xca 



Iv+k+l/2{x^/W^^c) 



^/J^^^c) ^\./JJ^J kW{v + k-l) 

When V = 0, a = c = ^jp then f{t) = ^-^—jj^r-----^ Io{^^) , so that if r ~ / then i/r has 
Rice distribution with parameters (a,/3). 

Example 6: Suppose = n and p„ = (^^)"~| '"(/j'')" ^^ /-^ where (x)„ is the Pochhammer 
symbol (see ( |46| ) in Appendix). Furthermore, we assume that the real valued sequences 
{aj}j=i,...,r and {/3j}j=i,...,g are such that pn > for all n > 1 and YIT Pn = K < oo. Then / 
and g are given by 

f{t) = — ^ rF,(ai, ...,ar;/3i, •••,/3g;ct) 

_ ce-^* ^ Kpn{ctY 



oo ^ X n-1/2 

,(x) = V2^e-- (t^J , 
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where rFq is the generahzed hypergeometric series (see Gradshteyn and Ryzhik (2000) 9.14) 



When r = q = 1 then p„ 



Mi 



— . Furthermore, ifa>0, (3 > a-\-l, then p„ > for all n and 



the series Yl^=iPn converge by Raabe's convergence test: 



hm n{pn/Pn+i 

njoo 



1) = lim n{(3 — a)/{a + n) = P — a > 1 . 

nfoo 



In this case / is given by 



fit) = ce-"' iF,{a,P;ct)/K, K = J^Pn 



where iFi is the confluent hypergeometric function of Kummer (see Gradshteyn and Ryzhik 
1(2000)1 9.21). 



We saw in Section 2 that for any strictly positive boundary, b{t) > 0,Vt > 0, with corresponding 



unconditional density fb{t), the two densities /o and fb are related as in (20) for all t. We will 



use this result to provide a connection between unconditional densities to the zero boundary and 



unconditional boundaries in the class C. To this end, let JC be the linear integral operator in (20) 
with b{t) = fit, i.e. 



JC :-- 



^{t - s)3/2 



Define the new class of densities := /CC, i.e. fo £ if there exists f G C (see (23)) such that 
/o = In this way, is a mapping which perturbs the probability measures of the class of 

densities C. Suppose Txb f ^ from Theorem 



4.1 



the matching density g is given in (24). Define 



the target density /o := ICf which is viewed as a target density for the zero boundary hitting time 



rx,o- By integrating equation (19) (with b{t) = fit) with respect to g, one finds that g is a solution to 



the RFPT for the zero boundary and target density /q. Therefore the pairs {fit, f{t)) and (0, /o(i)) 
result in the same matching density g. This result can succinctly be stated as follows. 

Corollary 4.1. Let f € C and suppose rx,o ^ fo '■= ^/ ^ C^- Then the matching density for the 



zero boundary and unconditional density /q is provided by (24)- 



Note that neither of the classes C and includes the other. For example when / is the density 
of the exponential distribution, there exists a unique g which matches the pair {fxt, /), for > 0, 
since the class of Gamma distributions includes the exponential distribution. However, in Section 2 
we saw that when the boundary is zero and the target density is exponential the unique solution to 
the RFPT Fredholm equation is g{x) = \/2Asin(\/2Ax) which is not a density function. 
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5. Boundary Transformations 

Suppose there exists a random variable which is the solution to the RFPT for the boundary h{t) 
and target density f{t). The question that we address in this section is the following. "How does the 
target distribution change as the boundary is modified, while the distribution of X is left unchanged?" 
We will put particularly emphasis on the linear boundary case. 

5.1. Afflne Transformations 

In this section we investigate affine modifications of the boundary. To this end, denote the 
boundary specific stopping time 

ni^t)-={t>^; X + Wt<b{t)] , (27) 

where X solves the RFPT for the pair {b{t) , f}ji^f^{t)) and has a density function g. Denoting the 
density of Tf,(() by f\,{t), from Ml) (conditional on X = x) we obtain: 



oo ;>oo 



'0 ^0 
Multiplying by g{x) and integrating over x, one finds the following connection between fi,(^t)+xt 
fbit) '■ 



e-"''(*)-*("'+"^/')Aw+At(t)rfi = K«) = / e-"^W-*"VV,(,)(t)dt . (28) 
'0 Jo 

This result relates a boundary specific transform of the density of Ti,^t)+\t to a boundary specific 

transform of the density of Tb(i) and the Laplace transform of the random starting point. For certain 

boundaries, the boundary specific transform is particular simple and we investigate the linear case 

next. 



Equation (28) reduces to a particularly simple expression when the boundary is linear. In 



particular, set b(t) = fit and apply an affine transformation to change the slope from fj, to i^, 



i.e. choose A = {i/ — fi), then (28) takes the form 



Utiau + aV2) = ~g{a) = Ut{afi + 0^/2) • (29) 



The above equality holds when a > max(— /i, — z^) so that equation (12) is valid for both linear 



boundaries fit and ut. Equation (29), therefore allows one to represent the Laplace transform of the 



density of the FPT to the new slope in terms of the density of the FPT of the old slope. 

Lemma 5.1. Suppose X solves the RFPT problem for the pair {fit, f^t{t))- Then, keeping the 
distribution of X unchanged, the stopping time of the randomized Brownian motion to the drift 
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adjusted boundary b{t) = vt has a density function f^t with a Laplace transform 

fut{s) = Ut - /i) + s + V2{fi -u)y^s + 1/2/2) , s>0. (30) 



Proof Let s = au + Then afi + = - fi) + s + \/2{fj, - u) ^/s + and using (|29| 

we obtain the result. □ 



Note that when s > the quantity in the brackets on the right side of ( 30 ) is positive for all // > 
and 1/ G M and thus the right hand side exists for any density function However, for /x < (30) 

is not a proper distribution. 



Interestingly, equation (30) allows one to determine the distribution of the FPT to the new 
boundary without knowledged of the initial density g of the starting point X. All that is required, 
is that g exists and it solves the RFPT problem for the boundary fit and the unconditional density 



ffj,t. In the specific case when the density /^^ is Gamma(a,/3), the relation (30) then reads 
Ut{s) = {1/ay (^^ + u{u-fi) + s + V2{fi - u)^s + uy2^ 

A + as + bVV+~c)'^ , (31) 



where A = I + au{v — /u), B = a\/2{fi — u), and C = 1/^/2. In Figure ??, we numerically invert the 
Laplace transform and plot f^t when /i = 1 for a few choices of u. This idea of changing drifts has 
potential applications in the context of finance and insurance, where such a drift change corresponds 
to changing probability measures from the historical measure to a valuation measure. 



As a final result for affine transformations, through Lemma A.l in Appendix A, we find the 
following explicit representation for f^t when jjL> u: 



This result may be more useful than the laplace transform representation (31) in certain cases. 

5.2. Scaling of Boundaries 

In this section we analyze transformations induced by a scaling of the boundary. In particular, 
we assume the boundary bx{t) is indexed by a parameter A > and satisfies the following scaling 
property 

bx{t) = bi{\h)/\ . (32) 

Note that the linear boundaries b\{t) = Xt satisfies this scaling property, further, scaling corresponds 
to modifying the slope. We are interested in how the stopping times 

Tx,x = inf{t >0;X + Wt< bx{t)}, X>0 , 
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indexed by the scale parameter, are related to one another. We will see that, through the scaling 
property of Brownian motions, we can investigate randomizing over the scale parameter rather than 
the starting point of the Brownian motion. Such a randomization can be viewed as a generalization 
of randomizing the slope of a linear boundary. In the RFPT problem we encountered the class C 



(see (23)) of mixtures of Gamma distributions. This class has a restriction on the minimum scale 
parameter in relation to the slope of the linear boundary. Thus, it is plausible that the randomization 
over the boundaries scale parameter may allow this restriction to be removed or at least modified. 
In the end, we hope to be able to probe a larger class of target densities. In light of these comments, 
we define the following modification to the RFPT problem. 

Definition 5.1. [Randomized Scaling First Passage Time Problem (RSFPT).) Given a scale class 
of boundaries h\ : [0, oo) — > M, and a probability measure ^ on [0, oo), find a random variable X such 
that /i is the law of the randomized scaling first passage time ti^\x- 

To solve this problem we assume, as usual, that the boundary is regular in the sense that P{tx,x = 
0) = for all values of A > and X. As a consequence of the scaling properties of Brownian motion 
we have the following easy Lemma. 

Lemma 5.2. If X has no probability mass at zero, then 

Tl,XX = J^Tx^X . (33) 

Proof. First, conditioning on X ^ 0, and using the scaling property of Brownian motions and the 
class of boundaries, we have 

^ = inf {t/X^ >0; X + Wt< bx{t)} 

= inf {t/X^ > ; 1 + Wt/X < bi{X^t)/{XX)} 
= inf {u > ; 1 + Wu< bi{X^X^u)/{\X)} 
= inf{'u>0 ; l + Wu<bxx{u)} 

Since the equality holds for every X ^ 0, and X has no mass at 0, we have the result. □ 

Suppose we have a solution gx{x) to the RFPT for the boundary /density pair {bx{t), f\{t)), i.e. 



Tx,x lias density fx when X has density gx- Then, through Lemma 5.2, the r.h.s. of (33) is known, 
and therefore the randomized scaling first passage time ti^xx is known and we therefore have a 
solution to the RSFPT. 
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Lemma 5.2 provides a solution to another related FPT problem where the intercept and the slope 
are both randomized simultaneously. In particular, consider the stopping time tx,^ (for a > 0) 
with X having no probability mass at zero - we can call the problem of finding the distribution of 
X given the distribution of ^x,^ the randomized double scaling FPT problem (RDSFPT). To solve 
it, consider the three random variables InX, Inrx,^ and Inri^^ and denote their mgf's by g, f and 
h, respectively. Then we have the following result. 

Corollary 5.1. If X has no probability mass at zero, then 

ff(2a) = £^ (34) 



Proof. From Lemma 



5.2 



we have Ti^a = x^'^^.f • Consequently, log(Tx,^) = 21og(X) +log(Ti,a); 
furthermore, log(X) is independent of log(Ti^„). This independence neatly implies ( |34[ ) for the mgf's. 
□ 



Corollary 5.1 implies that if we know the mgf of X then we have explicitly the mgf of tx,-^ and vice 
versa ~ assuming the mgf of ri^a is known. If those mgfs are invertible we have completely solved 
the RDSFPT. 

6. Conclusions 

In this article we introduced a problem related to both the forward and inverse first passage time 
problems. The new problem, coined the randomized first passage time problem (RFPT), seeks the 
density g of the initial starting point X of a Brownian motion such that the density of the FPT tx 



to a fixed boundary 6 is a given target function /. We prove two uniqueness results (see Prop. 2.1 



and 2.2) and provide an existence result (see Thm. [T]) for g. Furthermore, when the boundary is 



linear we provide several explicit examples of the density g based on the target / being a mixture 
of gamma distributions. Finally, we provide examples of how the target densities for boundaries in 
an affine and scaling class are related to one another. 

There are a number of theoretical directions left open. The first being making the existence 
result easier to check for boundaries other than the linear and square-root boundaries. The second 
is making the connection to Skorohod's embedding problem more explicit. We have already made 
progress on this direction and will report the results in another article. A third theoretical direction 



is to demonstrate the consistency of the existence results in Propsitions 2.1 and 2.2 

These results have several applications ranging from finance, to cosmology, to biology. In the 
financial context, the Brownian motion can represent the log of the leverage ratio of a company. 
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When this leverage ratio hits a critical barrier (the boundary), the company defaults. The dis- 
tribution of the default time is observable through yields of certain financial instruments (such 
as bonds and credit default swaps). Our results then provide a methodology for obtaining the 
distribution of initial log leverage ratio consistent with market prices. As well, the affine and scaling 
transformations allows one to change measure from statistical to risk-neutral (or pricing) measures 
in a natural manner. 
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Appendix A. Proofs of Results 



Proof of Lemma 2.1: Suppose h{t) > for some t > so that for y = —x the condition 
y = —X < b{t) — X is satisfied. Then the kernel of ([9| becomes free of x. Multiply Q by g{x) G G, 
assuming that there exists such a g which solves ([2]), and integrate x on (0, oo) to obtain: 

/■°° e-^^HJx/V2t) , , , /-"^ /■* e"^ / b(s) \ ,,,,,,, , 



Next we examine the right hand side of (35) and justify the exchange of the order of integration 



by the use of Fubini's theorem after we show the quantity under the double integral is absolutely 



integrable. Thus, using (43), we have: 



JQ 



b{s) 



[t-s)y 

6(a) 



j^^^^^^^f[s\x)\g(x)\dsdx (36) 



/■°° /"* (b(a))2 e'"V(*-=) 
< j^e ■^^^-^^j^—^^-^f{s\x)\g{x)\dsdx (37) 



1 -J^ n 



; (,_,)(n+l)/2 fis\xMx)\dsdx 



(38) 



where qn = ^2 [n/2j and z„ = 2"/2-K2J(n!/ [n/2j!) and 4 = 2"/2(n!/ [n/2\\). Let p„(s,t) = 

-i(-|iL.^,„)2 

(f-s)("+i)/2 — • have < Pn(0, t) < and Pnit, t) = since b is continuous. Thus Pn{-, s) 

is continuous and since it is finite at and at t then it is bounded on the interval [0, t] by, say, Mn{t). 
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Then, continuing from (38) we have: 



z'n 1^ _ s)("+i)/2 fi^\^)9{x)dsdx < z'^Mnit) f {s\x)\g{x)\dsdx 

/•oo 

= z'^Mn{t) / F{t\x)\g{x)\dx 
Jo 

where F{t\x) is the conditional cdf. Since b is continuous, for each t > 0, we can find an k{t) < 
and N{t) > —k{t) such that b{s) — x < —k{t) — x < for ah s < t and x > N{t). Therefore, for 



X > N{t), F{t\x) < 2^> 



x+k{t) 



< 



^g-(x + fc(t))2/(2t) 

x+k(t) 



and continuing from the last equality we obtain: 



z'^Mnit) / F{t\x)\g{x)\dx = 
Jo 

r-N{t) f-OO 

F{t\x)\g{x)\dx+ / F{t\x)\g{x)\dx 

lO JN{t) 



= z'nMnit) 
< Z^Mn{t) 



N(t) roo g-(x-+fc(t))V(2t) 

\g{x)\dx + Vt — — - — \gix)\dx 

JN{t) X + k(t) 



The integrals in the square brackets are both finite for all (7 G G (and more generally for all g ~ 
©(e^*^), < r < 2 for large x and which are absolutely integrable in the neighborhood of zero) and 



hence we can exchange the order of integration and the right side of (35 ) becomes a„(t). Furthermore 



the left side of ( 35 ) can be written as 

which is the Hermite transform of g{u^2t). Thus, since the Hermite functions form a complete 
orthogonal basis in the space L^(e~"^), if the Fredholm equation ([2| has a solution g ^ L^(e~"^), 



/•oo 

/ e~'^^ Hn{u)g{u\/2t)du , 
Jo 



it is unique (by the uniqueness of the Hermite transform) and it is given by (10) provided that h is 
continuous and b{t) > for some t. □ 



Proof of Lemma 3.1: Let 6 be a continuous function. Using (42) and (43) we have 

\£^eM-{o^Vm + m/V2t)^}\ = i(-7v^)"e-(-^)V„(M^) 

< i^„t"/2g-(^)^+2yK^^ 



where Kn = 2-L"/2j(„!/ [n/2j!). Thus 

^00 

Jo da^' 



-ab(t)-a^t/2 



\f{t)dt 



da^ 



f{t)dt 



< Kn 



^n/2^~ab(t)~aH/2+2,/[^\'-ii^ 



f{t)dt 
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Let e„ be such that t''/2^-abit)-aH/2+2,/[^\a^ ^ ^ ^ ^ rj.^^^ ^-^^ ^^^^ integral is 

bounded by 



Since b satisfies liuifiooiK^) + cut) > —oo then if Hmt|oo b{t) < 0, it follows that b{t) ~ —f^, 7 < 1 
for large t and the power in the exponent of the second integral above is dominated by — 
Similarly, if lim^ioo b{t) > then ab{t) dominates b{t) / \f2t. Therefore, the second integral is finite 
for all a > 0, n > 1, since 

lim i"/2g-"f'W-«^*/2+2VR^'^ = 

tfoo 



and b is continuous. Under the hypothesis of Lemma |3.1| it follows that 

,00 



'da" 



-ab{t)-c?tl2 



\f{t)dt < 00 



and therefore 



00 







n/2 



^ab{t)-aH/2TT ( Kt) + 



Complete monotonicity of r(a) is thus equivalent to the last integral being nonnegative for all 
a > 0, n > 1. This completes the proof. □ 



Proof of Corrolary 



3.1 



When b{t) = fit and f{t) = ^vw^^i fi,a,b > 0, the hypothesis of Lemma 



3.1 hold and thus we only need to check the complete monotonicity of r{a) by showing that (16) is 



nonnegative for all a > and n > 1. In this case (16) becomes 



_ ^ ["^ n+2h-l-x■^{a'+2a^l+2/a)/2TT ( [oi + fj)x 



dx 



The last quantity is strictly positive for n = 1 and all a > since Hi{z) = 2z. Thus, set n = 
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2k + 6, k > 1, where is or 1. Then, using Prudnikov (1986) (2.20.3(4)), we obtain: 



r(6)a^ 7o 
2 



2k+5+2b-l -x^(a^+2aii+2/a)/2 



H2k+5 



(a + fi)x 



dx 



r(6)of 



Cfc(a)(-l)^(l/2 -k- b)kT{k + 6 + b)x 



X2Fi{k + b + 6/2, k + b+{6+ l)/2; b + 1/2; 



- 2/a 
(/i + a)2 



Ck{a)r{k + 6 + b){-lf^{k + b- l/2){k + 6 - 1/2 - 1)... x 



r(6)af 



x(fc + 5-l/2-(fc-l))2Fi(fc + 6+^,fc + 6+^;fe+^; ^^^^/; 



where {.)k is the Pochhammer symbol (see ( |46| ) below), 2-^1 is the Hypergeometric function (see (45) 
below) and Ck{oi) = ^^^_^'^^2fc+2t+i • Thus, we see that (16) is nonnegative whenever 2F1 > 0. This is 



the case when ^'^ — 2/a > 0. Finally, the existence of X is guaranteed by Theorem [l] This completes 
the proof. □ 



Lemma A.l. Suppose v,d,B,C > are positive constants and A £ M and such that d + s + A + 
B\/ C + s > for any s > 0. Then the Laplace transform of 



h{t) := -p7^ l\t-xr^/'x^e-^^'-''^-''^'^^^-^dx 



2r{v)^l 

is given by 

h{s) = (d + A + s + BVs + C 
for s > 0. 

Proof We compute the Laplace transform directly: 

h{s) = r e-'%{t)dt= ^ r x-^e-^^-^^^'^ r u~^'\~<''+'^-'^ dudx 

^' Jo ^' '^nv)V^Jo Jo 

Tj foo / ^2Ty2 \ -1/4 

= / x^e--('^+^+^) K^y2 (xBVC + ^) dx 



B 2^/7^ 



00 

X 







2r(i;)v^ B 



D-lg-a::(d+A+s+BVC+s)^^ 



where we have used equation (3.471(9)) from Gradshteyn and Ryzhik (2000) | in the second line 
above. 
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Appendix B. Specicd Functions 

In this Appendix, we collect several useful formulae concerning the parabolic cylinder, Hermite 
and Hypergeometric functions. 

B.l. Parabolic Cylinder Function 

The parabolic cylinder function can be characterized as the solutions to the following differential 
equation indexed by p 

(Pu ( , 



^^, + ^p+l/2--J. = 0. (39) 

The solutions are u = Dp{z), Dp{—z), D-p-i{iz), D-p-.i{—iz). These functions also admit and 
integral representation for p < 0: 

Dp{z) = — — / e-'^'-^^'^x-P-Hx . (40) 
^ \—P) Jo 

As well, the parabolic cylinder functions for integer parameters are related to the Hermite polyno- 
mials 

Dn{z) = 2-^/'e-^y'H^(^-^^ . (41) 
Here iJ„ is the Hermite polynomial of degree n. 

B.2. Hermite Polynomials 

The Hermite polynomials are defined as the re-scaled n-fold derivative of the gaussian density as 
follows 

H^x) = i-ire^' {e-^') . (42) 
Some useful bounds on the Hermite polynomials are 



\Hnix)\ < 2"/2-L"/2J(„!/ Ln/2j!)e2^^L^ . (43) 

Finally, the Hermite polynomials form a complete basis on the Hilbert space L'^{R, e~^^) and satisfy 
the orthogonality relation 

Here 5„(m) = 1 if m = n and zero otherwise. 

B.3. Hypergeometric series 

The Hypergeometric function is defined as: 

,F,(a, 6; c; .) = 1 + 4x + '±±1}!^.' + ... = f ^jM^ , (45) 
c.l c(c+l).1.2 ^ [cnjnl 

n=0 
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where the Pochhammer symbol (z)„ denotes 

{z)n := z{z + l)...{z + n - 1), (n = 1, 2, ...), {z)o = 1 . (46) 
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